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Abstract 

The ambiguities and inconsistencies in previous thermodynamic treatments for the quark mass 
density-dependent model are addressed. A new treatment is suggested to obtain the self-consistent 
results. A new independent variable of effective mass is introduced to make the traditional ther- 
modynamic calculation with partial derivative still practicable. The contribution from physical 
vacuum has been discussed. We find that the properties of strange quark matter given by quark 
mass density-dependent model are nearly the same as those obtained by MIT bag model after 
considering the contribution of the physical vacuum. 
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I. INTRODUCTION 



It is generally accepted that the effective masses of particles will change with density 
due to medium effects. Many theoretical considerations, including the finite temperature 



quantum hadron dynamics (QHD) model|l|, Brown- Rho scaling[2j, finite temperature quark- 
meson coupling (QMC) model 3!], and etc., had been suggested to investigate the effective 
masses of mesons and nucleons. Besides theoretical study, many experimental results which 
predict the changes of particle masses with density have been shown. In particular, the 
experiments of TAGX collaboration had shown directly that when the density of the nucleon 
medium equals to 0.7no where no is the saturation density, the effective mass of neutral p- 
meson reduces to 610 MeV[4]. Both theoretical and experimental results confirm that the 
medium effects are important for studying the nuclear or quark systems. 

To illustrate the medium effects more transparently in theory, instead of first principle 
calculation, many authors introduced different hypothesis to represent the medium contri- 
butions, for example, supposed the density-dependent vacuum energy B(p) to modify QMC 
modelja, 1 suggested the density-dependent NNp coupling to address liquid-gas phase 
transition 7j, and etc.. Employing these hypothesis, many physical properties of nuclear 
matter, quark matter, nucleon system and hyperon system had been discussed. In particu- 
lar, to simplify the calculations, many ideal quasiparticle models in which the effective mass 
depends on the density and/or the temperature had been suggested for studying the quark 
gluon plasma (si, sj], gluon plasma [h| and strange quark matter 

The quark mass density- dependent (QMDD) model llj is one of such candidates. It was 
first suggested by Fowler, Raha and Weiner. In this model, the masses of u, d quarks and s 
quark are given by 

m*, = - — , (q = u,u,d,d), (1) 
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where B is a constant and m s o is the current mass of the strange quark; and ps, the baryon 
density, is defined as 
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where gi is the degeneracy, i (i = u, s, d) and i correspond to the quarks and antiquarks, 
respectively. The ansdtze in Eqs.(l) and (2) corresponds to a quark confinement mechanism 
because if quark goes to infinity, the volume of the system tends to infinity and the baryon 
number density goes to zero, then m* approaches to infinity according to Eqs.(l) and (2). 
The infinite mass prevents the quark from going to infinity 
Employing QMDD model and considering weak processes 



u + d <->• u + s, 



u + e + u e 



d — > u + e + u e 



u + e~ 



d + u e 



(5) 



and the condition of charge neutrality 



2Ap u = Ap d + Ap s + 3Ap e 



many 
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physical properties of strange quark matter [12J, [l3|, [14|, [l5|, [l6|, llll, UM, USA, 120, [2l|, [22, 
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25j and strange quark star [26|, [27| have been studied, and the results are in good 



agreement with those given by MIT bag model 28]. 

Although the density- dependent quark masses Eqs.(l) and (2) can mimic the quark con- 
finement mechanism, when we discuss the thermodynamic behaviors of the system with such 
quarks, many difficulties will emerge. The dispersion relation for quark 



e(k,p) = [fc 2 + m *(p) 2 ] 1/2 



(7) 



will make many usual thermodynamic relations with partial derivative no longer satisfied. 
As is well known in thermodynamics, a proper choice of independent variables will have a 
suitable characteristic thermodynamical function, from which all the thermodynamic quan- 
tities can be obtained by partial derivatives without integration. For example, with variables 
temperature T, volume V and chemical potential //, the characteristic function is the ther- 
modynamical potential Q = Q(T, V,/i). From the differential relation for reversible process 



dtt = -SdT - pdV - Ndfi, 



(8) 



we have 

S H4; p= -[nL' N= -w) T y (9) 

where S, p and iV are entropy, pressure and average particle number respectively. Other 
thermodynamic quantities such as internal energy U, Helmholtz free energy F, enthalpy H, 
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Gibbs function G, etc., can be calculated by the combination of the quantities we obtained, 
based on their definitions or relations. 

But for the QMDD model, Q is not only a function of T, V, p, but also depends ex- 
plicitly on the quark density p because of Eq.(7), Q = fl(T,V, p,m*(p)). How to tackle 



the thermodynamics self-consistently is stil 
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251 ]. The difficulty comes 



from the first and second laws of reversible process thermodynamics expressed by Eq.(8) 
and the partial derivatives by Eq.(9). Obviously, some extra terms involving the derivatives 
of m* will emerge when partial derivatives are calculated following Eq.(9). Unfortunately, 
these extra terms for different treatments in different references contradict each other. For 
example, for ideal quark gas system of quasiparticle with effective quark mass m* = m*(p), 
the pressure and energy density e were given by 
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in Ref.[18|], respectively. The ambiguity arises from the variable p in m*, because it is not 
one of the characteristic variables of thermodynamic potential. 

This paper involves from an attempt to clear above ambiguity and suggests a method 
to calculate the thermodynamic quantities from partial derivatives self-consistently. In next 
section, we will address the traditional thermodynamic treatment with partial derivative 
for QMDD model transparently and emphasize that the difficulty can not be overcome by 
usual method. In fact, the thermodynamic inconsistency for an ideal quasiparticle system 
with effective mass m*(p,T) to describe the quark gluon plasma or gluon plasma had been 
pointed out by many authors previously {§], O, Q . To avoid this difficulty, in Sec. Ill, we argue 



that we can calculate the thermodynamic quantities from equilibrium state and show that 
if we choose the quasiparticle effective mass m* as a new independent degree of freedom, 
we can calculate the thermodynamic quantities along reversible process by usual partial 
derivative self-consistently. In Sec. IV, we will apply this new treatment to study the QMDD 
model. We show that we can get rid of previous difficulties with partial derivative and take 
clear, reasonable and self-consistent results in the fresh, which coincide with the results in 
equilibrium state. We will consider the vacuum correction on the pressure and compare 
our result with that given by MIT bag model in Sec.V. We will show that even though our 
results are in agreement with that give n by Ref . 15|, [l6| , but the mechanism and the physical 



reasons are complete different. In Ref. 



15 



16] , the negative pressure comes from the partial 



derivative terms of m*. But in our treatment, it comes from the physical vacuum. Our 
mechanism is just the same as that of MIT bag model. The last section includes a summary 
and discussion. 
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II. INCONSISTENCY OF TRADITIONAL THERMODYNAMIC TREATMENTS 
WITH PARTIAL DERIVATIVE 



We first repeat the traditional treatment based on Eqs.(l) and (2). It has been established 



in Ref. 



15j, [16|, [13] systematically. In Ref.|17[], they derived 



(m\ _ ( d(vn) \ ( d(vn/N) \ _ ( d(n/ P ) \ _~ (m\ 

P ~ V 9V ) \d{V/N)) Tii - \d{l/p)) T ^ +P \dp) T ^ 

(21) 

Without losing generality, from here on we will display the formulae with only one component 
for simplicity and clarity. Eq.(21) is different from that of the normal expression p = 
— VL/V = —ft by the last additional term due to the density-dependence of effective mass. 
But Eq.(21) seems not so solid if we pay attention to the invariables. Noting that the 
chemical potential p is determined by the particle density constraint p = N/V, and the 
partial derivative of Eq.(21) is with respect to V but fixed T and p, or equivalently fixed T 
and p, in this process N must be changed. It is not a constant and so the third equal mark 
in Eq.(21) can not hold. The result in Eq.(21) is error. 

n 

To avoid this ambiguity, in Ref. [181], in stead of the thermodynamic potential f2, the 
authors introduced the Helmholtz free energy F to change the characteristic variable p to 
the average particle number N 

dF = -SdT - pdV + pdN. (22) 
By using F = Q + pN, they got 

After some calculations, they gave 

dQ (— dQ\ dp p dQ dm* 
P= ~dV~{ + l)p) dV + Vdm~*l)p~' { ' 

~ ~3T ~ \ + dp) dT dm* dT ' ^ ' 

For the internal energy, they required 



and found 



(28) 



But in their formulae, the invariable quantities for partial derivative have not been written 
down explicitly. To show the errors of their calculation, we write down these fixed quantities 
in the following by using the standard mathematical chain rule of partial derivative with 
composition function: 
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Eqs.(29) and (30) are just Eqs.(25) and (26) with the unchanged variable explicitly written 
out. Note that many partial derivatives are with four fixed variables. The term ( ) _ 

J K \ dm >' JV,N,n,T 

must vanish because m* = m*(T, p), p = N/V, to fix V, N and T is essentially to fix m*. It 
means all derivatives involving m* in Eqs.(29) and (30) should be zero, therefore the extra 

can not survive. In fact, according 
to the definition of Gibbs function G = pN — U — Ts + pV, we can directly find f2 = — pV. 
The extra terms involving derivative of m* will surely destroy this relation. 

The inconsistency dose not only appear in the expression of p and S, but also in N. If 
we calculate the particle number N from Eqs.(8) and (9), we obtain 
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where the summation is over all the quantum states with degeneracy gi. Since p = N/V is 
usually taken as a constraint to decide the value of p, ( 9m ^J'^ ) T y dose not vanish if m* 
depends on p explicitly, so the extra term in Eq.(31) modifies the particle number of the 
system. For an ideal quasiparticle system, the number of quasiparticle N = J2i9i n i, which 
is just the first term of Eq.(31), then we come to a conclusion that the result of Eq.(31) is 



incorrect. Besides, this result in Eq.(31) also conflicts with the requirement in Eq.(27), then 
the inconsistency of the traditional treatment with partial derivative is transparent. 

In fact, the inconsistency of usual thermodynamic treatment for ideal quasiparticle system 



has been shown by many authors 



25]. This section is not a new result. Our aim 



is to demonstrate the inconsistency for thermodynamic treatment with partial derivative 
along a reversible process explicitly. 

III. A THERMODYNAMIC ALLY CONSISTENT TREATMENT 



In Sec. II, we have shown the inconsistency of traditional thermodynamic treatment with 
partial derivative along reversible process. In this section, we will suggest a consistent ther- 
modynamic treatment for quasiparticle system with density and/or temperature dependent 
particle mass m*(p,T). We want to point out that our method is universal and QMDD 
model is just a special example for using this method. Our method can also be applied to 
discuss the thermodynamic behavior of the quark mass density and temperature dependent 
(QMDTD) model |l8 
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24| or other ideal quasiparticle system. 



Noticing that at a fixed instant of reversible process, the system is at an equilibrium 
state. Denote the temperature and density of this system as To and po, respectively, then the 
effective mass of the quasiparticle becomes constant m*(To,po) = m o- The system reduces 
to a usual ideal gas system with constant mass mo quasiparticles. For this equilibrium state, 
the corresponding thermodynamic quantities can be directly obtained. For example, for the 
case of one component Fermi system: 
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where rij is the particle number of the zth state and ^ is the corresponding degeneracy. 
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It can be clearly seen that these formulae are the same as those of the standard ideal 
gas. In fact, for example, if we derive the energy from the very beginning of the statistical 
definition in grand canonical distribution, we have 

U = E = I^r^e"^-^ 
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where the invariables of \i is clearly seen, while m* and V are fixed since En is a function of 
m* and V for each N. Take the quasiparticle ideal gas thermodynamic potential, Eq.(33), 
which means the grant partition function satisfies 



InS = $>ln(l + e -/V™* 2 + fc ?-M)), (39) 



it is easy to find 
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where rij = e/3[e(m<<j |,. ) _ M]+1 is the average particle number of the ith state. This result consists 
with the interaction-free quasiparticle picture of the QMDD model. 

We see, from Eqs.(32)-(37), the contribution of medium effect is included in the effective 
value of mass, and appears in the exponential of the Fermi distribution. A remarkable 
property of these formulae is that the extra terms related to the partial derivative of m* 
do not appear. This is reasonable because these thermodynamic quantities are functions of 
equilibrium state, they do not depend on the change of the quasiparticle mass, since at this 
state the quasiparticle mass m*(T,p) becomes constant mo(T Q ,p ). 

To compare above treatment with those given in Sec. II, as an example, we calculate the 
entropy of the ideal quasiparticle system. Denote the entropy calculated by equilibrium 
state as S sta . From Eq.(34)-(37), we have: 
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(41) 



This is a familiar formula for equilibrium state whose physical meaning is transparent. While 
we denote the entropy calculated by the partial derivative of Q following Eqs.(8) and (9) as 
Sder, we have 
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and for the equilibrium state with m*(T ,p Q ) = m , 
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(43) 



The first term of Eq.(42) is just the result given by equilibrium state. The difference 
between these two treatments is significant. They can not be accorded together. Noticing 
that the contribution of medium effect at equilibrium state is included within the value of 
m* in the distribution, and the entropy describing disorder of quasiparticles in a system dose 
not depend on the intrinsic quantity such as effective mass of quasiparticle, the correctness 
of S s ta is obvious. The thermodynamic consistency of the new treatment is transparent 
because it is based on an equilibrium state and m (T ,p Q ) is constant. 

If one hopes to extend this treatment to a reversible process, the accordance between the 
results from calculation along a reversible process and results obtained at a fixed equilibrium 
state must hold. So in order to get a self-consistent calculation between equilibrium state and 
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reversible process in thermodynamics, we are compelled by above discussions to introduce 
an intrinsic degree of freedom m* for quasiparticle in thermodynamic system to describe 
medium effect. We should rewrite Eq.(8) as: 



dtt = -SdT -pdV - Ndp + Xdm\ (44) 



then Eq.(9) becomes 



\ dT ) V ,^ ' \ dV )T^ \ d V) T ,V,m*' ' \ dm * J T,V^ 

(45) 

where X is an extensive quantity corresponding to the intensive variable to*. In Eq.(44), the 
intrinsic degree of freedom m* for quasiparticle has been added as an independent variable in 
thermodynamic system. The quantities S, p and N shown in Eq.(45) are in agreement with 
the results obtained by the formulae of equilibrium state in Eqs.(32)-(37), because we have 
fixed to* as an unchanged parameter in the partial derivative calculations. For example, the 
pressure given by Eq.(9) is 



( dn(T,V,fi,m*(p B )) \ 
P ~ { dV ) 
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while given by Eq.(45) as 



/ 7> 



(46) 



p= ~{ — ov — T r~v- (47) 

In Eq.(46), the first term is just the result of equilibrium state, but the second term, ac- 
cording to the treatment in Sec. II, will lead to inconsistency in thermodynamics. 

Usually the thermodynamic parameters such as S, p, T ... depend on the whole system. 
They are independent with the intrinsic property of the particle or subsystem, no matter 
the subsystem is a simple point particle or a quasiparticle with inner structure and dif- 
ferent intrinsic properties. Ordinary thermodynamic variables depend on the collection of 
the subsystem only. Similarly, the mass is an intrinsic quantity of a particle, it does not 
affect on collective thermodynamic properties of the whole system. When we consider the 
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medium effect or the confinement mechanism, and summarize this effect into the effective 
mass m*(T,p) under quasiparticle approximation, the dynamic interaction can be concen- 
trated on the effective mass m* of quasiparticle by using the finite temperature quantum 
field calculation or directly by confinement anstdze. But the macro thermodynamic 
variables cannot describe these micro dynamic interactions. We must choose new variables 
to represent these dynamic interactions or the medium effect. Obviously, the effective mass 
m* appears as the suitable independent variable. Further more, we want to emphasize that 
the application of this new independent variable m* to rewrite Eqs. (8) and (9) is limited to 
quasiparticle approximation only. In quasiparticle approximation, the interaction between 
particles, the medium effect or the confinement mechanism are summarized in the effective 
mass m*. Then m* becomes an independent variable to represent all those physical effects, 
its dependence on density and temperature is due to equilibrium condition of the whole 
system. Introducing m* in quasiparticle physical picture to represent the medium effect and 
taking it as a variable is a twin in thermodynamics of quasiparticle system. If we employ 
finite temperature quantum field theory to calculate the thermodynamic potential beyond 



quasiparticle approximation 



31 



321 ] . we can use Eq.(2) to calculate the thermodynamic 



quantities and investigate the interaction between particles by different orders of Feynman 
diagrams. Then we will not need to introduce the quasiparticle and use the corresponding 
Eqs. (44) and (45). 



IV. QMDD MODEL 

To illustrate the basic difference between our treatment based on Eqs. (44) and (45) and 
those of traditional treatment based on Eqs. (8) and (9), we employ QMDD model as an 
example to calculate the thermodynamic property of strange quark matter. We fix the 



parameters in Eqs.(l) and (2) B = 170 MeVfm-^m s0 = 150 MeV, as in Ref.[2l|. The 
temperature is set at T = 0, as that of Ref. BBQ, for convenience of comparison. Our 
results are shown in Fig.l and Fig. 2. In Fig.l, we draw the curves of energy per baryon 
e/pb vs. baryon number density ps, where the solid curve refers to our treatment and the 
dashed curve see that the saturation point of the solid 

curve is ej ps = 906.3 MeV at ps = 0.433 fm~ 3 , which locates in a reasonable range. The 
saturation point of the dashed curve is ej ps = 1042.8 MeV at ps = 0.692 fm~ 3 . It is higher 
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than our result because it has an additional term on ( -J^- 

The equation of state for different treatments are shown in Fig. 2, where the solid curve 



refers to our treatment, and the dotted curve, dashed curve refers to that given by Ref.[l£ 
and Ref.|l7(], respectively. We see from Fig. 2 that there are two remarkable differences 
between the solid curve and others. Firstly, the pressure is always positive in our treatment. 
This is of course reasonable because our Hamiltonian H = Y.k,s £ ksO'i s a ks is an ideal gas 
Hamiltonian, it cannot give negative pressure in principle. But for the dashed curve and 
dotted curve, the pressure becomes negative in the small energy density regions due to the 
incorrect modification term 



M Hi - mf - m* In 



VTL: 
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PPb)t^} Y48tt 2 p B 

Secondly, although the tendencies of all the curves are similar in large energy density region, 
in small energy density region, their behavior are very different. In particular, the end 
points at ps — > for these three curves are different. They are {e = 170MeV fm _3 ,p = 
0} for solid curve, {e = 340MeVfm" 3 ,p = -170MeV fm" 3 } for dashed curve and {e = 
170MeV fm" 3 ,j9 = -170MeV fm" 3 } for dotted curve, respectively. 

Can we get negative pressure pressure in QMDD model, as that of the MIT bag model? 
We will answer this question in next section. 



V. THE CONTRIBUTION OF PHYSICAL VACUUM 

To answer this question, let us recall that in the MIT bag model the negative pressure 
comes from the physical vacuum. If we extend the quasiparticle Hamiltonian to 

H = £ k s at s a ks + H , (49) 

k,s 

where Hq is the system energy in the absence of quasiparticle excitations. Noting that in a 
quasiparticle system, the physical vacuum energy H is possibly a function of T and p since 
Eks depends on T and p via effective mass m*. Corresponding to Ho, the additional term of 
thermodynamic potential is Qq, and the pressure becomes 

P = y — =p-B , (50) 
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the energy density becomes e + B , but the entropy density s = (e + p — J2iPiPi)/T is 
unchanged. 

Now we turn to determine f^o i n QMDD model. By means of the vacuum correction to 
guarantee the thermodynamic consistency of quasiparticle system, it was first suggested by 



Gorenstein and Yang 
satisfy 



101 ]. they gave a set of constraints. While in Ref. 25J, Qq is chosen to 



o. 



V 9 PB / {M 

where fl = flo/V is virtually B in Eq.(50). At zero temperature, Eq.(51) reduces to 
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because m* depends on p B only, Cl can be obtained by integration. Based on Eqs.(l), (2) 
and (33), we obtain 
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Chosen (Woo) = £ as that of Ref. 
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fi (p 



B 



PB 



E 



QjBm* 
Wp% 



PiVPi ~ m 



f-mf\n 



Pi + ^P 2 i-mf 
m* 



dp B + B^. (55) 



Eqs.(53)-(55) are just Eqs.(7),(8) and (14) in Ref. [251] . then the stability and the negative 
pressure have been shown there. The negative pressure and the stable point dose exist, we 
do not recalculate the numerical results here again. 

Of course, the constraint Eq.(51) is one of the possible ways to get f2o only. The aim 
of this choice is to demonstrate that we can get negative pressure and a stable strange 
quark matter. We must emphasize here again that the negative pressure as well as the 
stability of the strange quark matter come from physical vacuum. It is completely different 



from that of Ref. 



•angequarK 



in which the negative pressure comes from the extra partial 
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derivative term ps-§^- As has been shown in Sec. II, this term cannot exist since it destroy 
the thermodynamic consistency Our mechanism is the same as that of MIT bag model. 



The conclusion [15J that the properties of strange quark matter in QMDD model are nearly 
the same as those obtained in the MIT bag model is still valid, but their argument was 
incorrect. The correct argument is shown by Sec.III-V. 

VI. SUMMARY AND DISCUSSION 

In summary, we have shown the shortcomings of the previous treatments for QMDD 
model, which are based on the partial derivative of thermodynamical functions along re- 
versible process. A new method is suggested. For the reversible process, we have introduced 
a new intrinsic degree of freedom m* for the system of quasiparticle. We have proved that 
the thermodynamic quantities calculated by the partial derivatives concerning this indepen- 
dent variable are in agreement with those obtained by the equilibrium state. The difficulties 
and controversies in previous references are removed. We also find that the properties of 
strange quark matter in QMDD model are nearly the same as those obtained in MIT bag 



model 23j], if we take the vacuum contribution into account. Our method is applicable 
if the quasiparticle with effective mass is introduce to represent the medium effect, since 
such quasiparticle models are quite commonly used in nuclear matter, quark matter and 
quark-gluon plasma, we hope our discussion might help the theoretical study in these fields. 
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FIG. 1: Energy per baryon as a function of the baryon density pb at T = for different treatments 
of the QMDD model. The solid curve represents our treatment and the dashed line refers to the 
treatment in Ref.[15||. 
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FIG. 2: Pressure P as a function of the energy density e = U/V for different treatments of the 
QMDD model. The solid curve, dotted curve and dashed curve represent our treatment, the 
15 ] and the treatment in Ref. \v\ . respectively. The tendencies of these curves 



treatment in Ref. 



are similar at large energy density region, but at small energy density region, different treatments 
have quite different behaviors. In our treatment, the pressure never goes to negative. 



18 



